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ABSTRACT 

The direct detection of dark matter on Earth depends crucially on its density and 
its velocity distribution on a milliparsec scale. Conventional N-body simulations are 
unable to access this scale, making the develo pment of oth e r app roaches necessary. In 
this paper, we apply the method developed in lFantin et al.l (|2008l ) to a cosmologically- 
based merger tree, transforming it into a useful instrument to reproduce and analyse 
the merger history of a Milky Way-like system. The aim of the model is to investigate 
the implications of any ultrafine structure for the current and next generation of 
directional dark matter detectors. We find that the velocity distribution of a Milky 
Way-like Galaxy is almost smooth, due to the overlap of many streams of particles 
generated by multiple mergers. Only the merger of a 10 10 Mq analyse can generate 
significant features in the ultra-local velocity distribution, detectable at the resolution 
attainable by current experiments. 

Key words: methods: numerical - Galaxy: evolution - Galaxy: halo - Galaxy: kine- 
matics and dynamics - Galaxy: solar neighbourhood - dark matter 



1 INTRODUCTION 

In the last forty years, dark matter (hereafter DM) has been 
one of the more challenging topics in astronomy. After the 
first pioneering studies (Zwicky 1933 and Smith 1936), con- 
vincing evidence for the existence of DM came from the ro- 
tation curves of spiral galaxies (Freeman 1970, Einasto 1974, 
Rubin & Ford 1970, Rubin et al. 1978, 1980, 1982, 1985). Its 
presence has been also supported by recent indirected mea- 
surements of its density, obtained by surveys such as SDSS 
(Tegmark et al. 2006), or by satellites like WMAP (Komatsu 
et al. 2009, Larson et al. 2010). Despite all this evidence, the 
nature of DM is still unknown. The most widely accepted 
idea is that it consists of non-relativistic, weakly interacting, 
non-baryonic particles (Cowsik & McClelland 1973; Szalay 
& Marx 1976), emitting no (or very little) electromagnetic 
radiation. 

Numerical simulations provide a reliable method of cal- 
culating the DM distribution on large scales, making robust 
predictions for its clustering. The resolution of the first sim- 
ulations (Peebles 1982, Frenk et al. 1985) was not sufficient 
to resolve the inner parts of the DM halos, but the gap has 
started to be filled in the last decade, with N-body simu- 
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lations able to provide results for halos of the size of the 
largest galaxy clusters (Springel et al. 2005, Boylan-Kolchin 
et al. 2009). More recently simulations such as Via Lactea 
II (Diemand et al. 2008), GHALO (Stadel et al. 2009) and 
Aquarius (Springel et al. 2008) were able to reproduce the 
formation and the structure of a Milky Way-like DM halo 
down to a resolution of about 100 pc. 

Particle physics provides us with various well-motivated 
dark matter candidates, including weakly interacting mas- 
sive particles (hereafter WIMPs). They interact with or- 
dinary matter through elastic scattering on atomic nuclei 
(Goodman & Witten 1985), and many experiments are cur- 
rently underway aiming to detect this phenomenon, such 
as CDMS II (Ahmed et al. 2010), XenonlOO (Aprile et al. 
2010) and Zeplin III (Akimov et al. 2010). The goal of these 
experiments is to measure the number of recoil events per 
unit energy and, in some case, its temporal or angular de- 
pendence. These quantities depend on the local DM density 
and speed distribution of the incident particles in the detec- 
tor rest frame (Jungman et al. 1996, Lewin & Smith 1996). 
Therefore, the possible presence of a significant amount of 
fine-grained structure in the DM distribution in the Solar 
neighbourhood, such as overdensities and streams, would 
have important consequences for these direct detection ex- 
periments. 
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The establishment of a reliable model for the Milky Way 
still represents a significant source of uncertainty. One of 
the main reasons is that the scale relevant for direct de- 
tection experiments (~ O.lmpc) is well beyond the resolu- 
tion of current simulations. Furthermore, DM particles are 
"cold", allowing the formation of micro-halos with masses 
down to 1O _6 M (Green et al. 2004). Unfortunately, such 
micro-halos are far smaller than the smallest subhalos re- 
solvable by simulations, which have mass ~ 10 4 M Q and size 
~ lOOpc. Consequently, other approaches are required to 
resolve the ultra-small scales probed by direct detection ex- 
periments. 

Stiff and Widrow (hereafter SW) calculated the DM ultra- 
local velocity distribution f(v) at a single spatial point of 
phase-space, identified as an ideal terrestrial detector (Stiff 
& Widrow 2003). This result was achieved using a reverse 
technique that allowed them to reach a resolution that can- 
not be obtained by N-body simulations. They found that 
f(v) is characterised by the presence of a small number of 
discrete streams of particles. Unfortunately, the reverse inte- 
gration is unstable, requiring the introduction of a softening 
length of 20 kpc into the adopted gravitational force law. 
This softening is large enough to affect the DM phase-space 
distribution imprinting on a terrestrial detector. 
Other theoretical investigations have argued that the ultra- 
local WIMP distribution consists of a much larger number 
of streams, of the order of 10 5 " 6 (Helmi et al. 2002; Vogels- 
berger et al. 2009). Overlapping DM streams in the inner 
part of the halo lead to more complex configurations, which 
may ultimately be indistinguishable from a smooth distri- 
bution. Moreover, direct cosmological simulations show that 
the most massive local streams only contribute about 1 per- 
cent of the local DM density (Diemand et al. 2008, Vogels- 
berger et al. 2011). Schneider et al. (2010) found that at 
Earth's location the DM ultra-fine distribution is composed 
of over 10 3 streams, which are generated by the tidal disrup- 
tion of DM "cold" microhalos (Schneider et al. 2010). They 
conclude that the material these streams are composed of is 
not dense enough to be relevant for detection experiments. 
Recently Vogelsberger and White (2011; hereafter VW) have 
argued that the picture on ultra-local scales is not too dif- 
ferent from that suggested by N-body simulations. Using a 
new technique for calculating the phase-space distribution 
function in the neighbourhood of a simulation particle (Vo- 
gelsberger et al. 2008), they estimate the presence of about 
~ 10 14 unresolved streams at the Solar position. Their con- 
clusion is that the local velocity distribution is smooth (Vo- 
gelsberger & White 2011). 

In this paper, we present a new method of simulating 
the DM distribution of the Milky Way, building on the ap- 
proach developed by Fantin et al. (2008; hereafter Paper I). 
Applied to a merger tree describing the history of a Milky 
Way-like galaxy, the model attributes to each progenitor of 
the halo a characteristic DM distribution, using the method 
developed in Paper I. The final result, obtained by summing 
the contributions of all the progenitors, provides a model of 
the ultra-fine DM distribution in the Solar neighbourhood 
at arbitrarily high spatial resolution, obtained without the 
computational overhead of a complete numerical integra- 
tion. 

The paper is organised as follows. In Section 2 we 
present the technique to calculate the ultra-fine DM dis- 
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Figure 1. The rotation curve of the isochrone potential, with the 
scale length taken to be equal to the Solar radius. Although not 
intended as a fully realistic model, it is a reasonable approxima- 
tion to the measured rotation curve of the Milky Way (Fig. 5 of 
Sofue et al. 2009). 

tribution of a merger tree. Section 3 contains our results, 
Section 4 draws some conclusions, and Section 5 provides a 
summary. 



2 THE MODEL 

In Paper I we modelled the interaction between an unbound 
system of particles and a Milky Way-like galaxy. Although 
not intended as a quantitative description of the Milky Way, 
this approach makes possible a detailed exploration of phase- 
space, and a thorough analysis of the expected signature of 
a series of merger events in a terrestrial DM detector. The 
method has the great benefit of being analytically soluble, 
allowing the calculation of the dynamics of the merger re- 
markably simply. Furthermore, the gravitational force does 
not have to be artificially softened, providing a significant 
improvement in accuracy. To describe the Milky Way grav- 
itational well we adopt the isochrone potential. Despite not 
being a realistic representation of a complex system such 
as the Milky Way, for example causing an overestimation 
of the granularity of the system at the present time, the 
qualitative properties of the potential are similar to those of 
our galaxy. Fig [T] presents the rotation curve calculated with 
a scale length equal to the solar radius. Both the value of 
the circular speed at the Solar radius and the overall shape 
of the rotation curve are in reasonable agreement with the 
Milky Way rotation curve shown in Fig. 1 and 5 of Sofue et 
al. (2009) 

Concerning the merging halo, we model its initial DM 
phase-space distribution by assuming a bivariate Gaussian, 

/(r,v) oc e -K*-*~) a /a«J)] e -[(v-v sat ) 2 /2 CT ?] ; (1) 

where (r sat , v sat ) is the initial location of the merging halo 
in phase-space, which allows us to simply parametrise its 
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main properties, where a s models its initial spatial extent 
and <7 V its velocity dispersion. 

Galaxy formation models (Kauffmann et al. 1993; Cole 
et al. 1994) and N-body simulations (Springel et al. 2005, 
Boylan-Kolchin et al. 2009) both suggest that galaxies 
formed though a continuous and regular hierarchical pro- 
cess. We therefore embed the model for a single merger into 
a cosmological context by combining it with a merger tree. 
This approach allows us to describe a system in which lots 
of DM satellites fall into the Galactic potential well. Fur- 
thermore, it produces a more realistic estimate of the final 
velocity distribution of the Galactic halo than the one pre- 
sented in Paper I because the merger tree describes a halo of 
1O 12 M0, which has not undergone any recent major merger, 
as in the Milky Way. The tree has been kindly generated by 
Andrew Benson, using a semi-analytic method (Cole et al. 
2000, Parkinson et al. 2008) based on the assumption of the 
following cosmological parameters: fio = 0.25, Ao = 0.75, 
fib = 0.045, ho — 0.73 and as = 0.9. The information pro- 
vided by the merger tree are the virial mass M of the merg- 
ing subhalos and the scale factor a at which they fall into the 
main progenitor. The mass resolution of the tree is 10 8 Mq. 
Knowing the epoch in the past at which each satellite falls 
into the Milky Way's halo and the initial position of the 
merging satellite in phase-space at that time, (r sat , v sat ), the 
model evolves the system analytically backwards in time. 
This allows us to map out the full velocity distribution 
within the "ideal" detector, which is located at the Solar 
position, at a distance ro = (8.5, 0, 0) kpc from the Galac- 
tic centre. We assume that each subhalo falls into the host 
system on a radial orbit from a random position on a sphere 
of radius r sa t = ?Vir, where r v i r is the (redshift dependent) 
virial radius of the Milky Way halo. Experiments reveal that 
these assumptions are not critical. For example, replacing 
the radial orbits with orbits that have a peri-to-apo-centre 
ratio of 1:6, the median found in cosmological simulations 
(Ghigna et al. 1998, Diemand et al. 2008), results in veloc- 
ity distributions that are indistinguishable. Similarly, start- 
ing the disintegration of each subhalo at its first pericentre 
passage rather than at the Milky Way's virial radius on its 
initial infall does not substantively alter the final velocity 
distribution. The quantity v sat is approximated by the ve- 
locity of a body falling from an infinite distance, 
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2GM 



Finally, we use the virial theorem 



Mvir = 



G 



and the definition of virial mass 

Mvir = -TrAvirPcritTvir ) 



(2) 



(3) 



(4) 



to relate the velocity dispersion, the spatial extent and the 
virial mass. The virial overdensity A v i r is defined as the den- 
sity relative to the mean density within r v i r times the critical 
density p cr u at that redshift (Bryan & Norman 1998), and 
Pcrit is the critical density. Assuming a value for the con- 
centration parameter c = r v i T /a s — 10 (Bullock et al. 2001; 
Benson 2005), Eqs. (0) - (|4]) give an estimate of the two 
characteristic scales of the merging subhalo. 
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Figure 2. Energy spectrum of a Milky Way-like halo for an evo- 
lution time of ~ 13.2 Gyr. The merger history of this halo of 
10 12 Mq does not include any recent major merger. One velocity- 
unit corresponds to 570 km s — 1 . 



To summarise, the principal steps of this modelling process 
are: 

• Pick a merger tree that gives a realistic representation 
of a Milky Way-like system. 

• Calculate the initial conditions for each subhalo in the 
merger tree. 

• Evaluate the velocity distribution, observable within 
the detector for each subhalo, using the model developed 
in Paper I. 

• Evaluate the total velocity distribution by summing all 
the constituent mergers. 



3 RESULTS 

The model we have developed above is not intended to pre- 
dict in quantitative detail the experimental signal that a 
terrestrial detector would observe. Nevertheless, it is inter- 
esting to define quantities related to the ones measured by 
directional detectors, with the aim of investigating qualita- 
tively the properties of the system. In Section 3.1 we first 
determine the local DM energy spectrum of a Milky Way- 
like halo as a quantity that can be directly measured by 
experiments. Secondly, we present the velocity distribution 
as a function of direction and speed, quantities that can 
give important indications about the evolution history of 
the Galactic halo. In Section 3.2 we present the contribu- 
tion to the velocity distribution produced by the addition of 
a new progenitor to the same Milky Way-like halo, in order 
to determine if late-merging subhalos might be detectable. 



3.1 Formation of a Milky Way-like dark matter 
halo 

The energy spectrum is an observable quantity, and the anal- 
ysis of its features can shed light on the evolution history of 
a system. For example, the particles of a stream produced by 
the recent merger of a subhalo normally group around par- 
ticular values of the energy, forming "overdensities" . Fig. [2] 
presents the energy spectrum of the Galactic halo for an 
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Figure 3. Fractional departure, r/, of the velocity-squared of the 
Galactic halo simulated in Figure [4] from a smooth distribution 
for an evolution time t ~ 13.2 Gyr. As in Fig. [2] one velocity-unit 
corresponds to 570 km s _1 . 



evolution time of ~ 13.2 Gyr. In the plot we see that the en- 
ergy spectrum does not show any small-scale feature. This 
smoothness arises because the halo is composed of a number 
of particles large enough to wash out any possible "overden- 
sity" produced by single streams. To determine the level of 
residuals from a completely smooth distribution, we com- 
pare the current configuration to that obtained by evolving 
the system to many times the current age of the Milky Way, 
thus ensuring that it is fully relaxed. The fractional depar- 
ture from this relaxed configuration is then given by: 

No - iVrol 



I] 



(5) 



where N e is the number of particles in a range of velocities 
at ij and N le \ is the number of objects in the same range 
once the system is relaxed. 

In Figure[3]we present r\ for the same system as Fig. [2] From 
the plot it is evident that the fractional departure is abso- 
lutely indiscernible from the distribution that the halo as- 
sumes when completely relaxed. The absence of any feature 
leads to the conclusion that the ultra-fine DM distribution 
in the Solar neighbourhood is composed of a huge number of 
streams. These overdense structures, clearly visible when we 
analyse the interaction of a single satellite, overlap, generat- 
ing a smooth DM distribution and consequently a featureless 
fractional departure. 

In principle, more information might be preserved if we 
look at the full velocity distribution as a function of angle, as 
recorded by some DM detectors. Fig. [4] is a representation of 
the velocity distribution /(r, v) of the Milky Way-like halo 
simulated above in the (cos 0,v) space, where 6 is defined 
as the angle at which DM particles enter the detector, mea- 
sured relative to the direction of Solar motion in the Milky 
Way. The plot is grey-scaled, and the shade of each pixel 
covers a range of six orders of magnitude in /(r, v). This 
range is large enough to show the most significant contribu- 
tions to the velocity distribution of the system. Assuming 
«max = 570 km s , a value slightly larger than the escape 
speed from the Milky Way (y esc = 544 km s -1 , Smith et al. 
2007), and w m i„ = km s - , the velocity resolution of each 
cell along the :r-axis is ~ 3.8 km s _1 . The angular resolution 



varies between ~ 10° (for small angles) and ~ 1° (for large 
angles). To investigate the presence of overdensities in the 
angle distribution we apply a boxcar smoothing technique. 
This removes all the small differences in the velocity-space 
distribution, mainly due to pixelation artifacts. Fig.|4]shows 
that in a system resembling our galaxy no clear features 
are present, except for a central overdensity. The overdensi- 
ties are relative to the surrounding pixels. The fact that the 
overdensities are not easily distinguishable means that a full 
quantitative analysis is not justified, but there is clearly in- 
formation in these structures. The particles composing this 
region have velocities in the range 285 km s _1 < v < 400 
km s _1 , which in the plot correspond to 0.5 < v < 0.7, and 
angles 60° < 6 < 120° (-0.25 < cos6> < 0.25). A second 
level of overdensities seems to be present in the same re- 
gion, though it is difficult to separate them from the main 
one. These weak features are due to recent and incomplete 
merger activity, and their existence raises the possibility of 
disentangling at least the recent merger history of the halo. 
This could be achieved by directional detectors with angle 
and velocity resolution of ~ 1° and ~ 3.8 km s _1 . 
An investigation on the origin of these two levels of over- 
densities can be carried out by analysing the contributions 
to the final velocity distribution of progenitors with similar 
mass. Using the same smoothing technique applied above, 
in Figure 5 we group together the satellites with mass of 
the same order of magnitude. The three groups we adopt 
are: 10 8 - 10 9 M (Fig. 5a), 10 9 - 10 10 M© (Fig. 5b) and 
10 10 - 10 11 M (Fig. 5c). The three figures show that the 
contributions from these groups of subhalos are heteroge- 
neous. The number of direct progenitors in the range of mass 
10 8 — 10 9 M Q is large (~ 10 3 ) and their total contribution to 
the final velocity-space distribution is not negligible. Most 
of the particles have 75° < 9 < 105° (-0.25 < costf < 0.25), 
but there is no evidence for a preferential range of velocity. 
This is due to the fact that these progenitors are compact 
and the timescale on which they get disrupted is longer than 
that of more massive systems. As a consequence, even after 
an evolution time equal to the Milky Way's age they are not 
completely disrupted. 

Looking at Fig. 5b and 5c we note that the contributions 
produced by subhalos of mass 10 9 " 10 M and 10 10 " 11 M 
are very similar. A significant difference can be seen: the ve- 
locity distribution produced by 10 9 Mq subhalos does not 
show any particular feature, while that from 10 10 M© satel- 
lites has a central overdensity. This suggests that the most 
pronounced overdensity of Fig. 3] is mostly produced by 



satellite with masses in the range 10 



M, 



3.2 Contribution of single merger events to the 
ultra-fine distribution of a Milky Way-like 
system 

This motivates the study of the contribution from single 
mergers. The method we adopt to investigate this issue con- 
sists of adding a DM satellite to the Galactic halo and of 
analysing its contribution to the final velocity-space dis- 
tribution. The properties of the additional progenitor that 
we vary are the mass and the time at which it falls into 
the Galaxy. For the mass, we consider four different cases: 
1O 7 M , 10 s Mq, 1O 9 M and 10 10 Af Q . The main reason 
that motivates the limit of 10 10 Mq is the absence of major 
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Figure 4. The smoothed velocity distribution of a Milky Way-like DM halo of mass 1O 12 M0 as a function of cos0 and v. The angle 
distribution is smoothed over a (5 X 5) box. The detector is located at a position ro = (8.5, 0, 0) kpc. We use an array of dimensions 
n v ,ng = (150 X 120). The dimension of a pixel on the x-axis is ~ 3.8 km s _1 . The angular resolution on the y-axis varies between 10° 30' 
(for small angles) and ~ 1° (for large angles). The greyscale, with the shade of each pixel determined by the corresponding array of the 
distribution function /(r, v), covers six orders of magnitude. Black corresponds to the upper limit and light grey to the lower one. The 
black circle highlights the region where overdensities, relative to the surrounding pixels, are present. 



mergers in the recent history of the Milky Way, based on 
theoretical considerations which show that the accretion of 
a 1 : 10 satellite galaxy could destroy the disk of the Milky 
Way (Purcell et al. 2009) . To get a complete overview of the 
evolution of the mergers we select four times, looking back 
into the past, at which the progenitor fell into the host halo: 

• ti ~ 1.75 Gyr ago (z = 0.14): the subhalo is making its 
first passage through the detector today at redshift z = 0. 

• t2 — 2.0 Gyr ago (z = 0.17): stage at which the satellite 
is not passing through the detector today. 

• £3 ~ 7.5 Gyr ago (z = 1.0): the merger is at an interme- 
diate stage of its evolution. 

• £4 ~ 13 Gyr ago (z = 8.5): the satellite merged at an 
early stage of the formation of the Milky Way. 

The results of these simulations are summarised in Table [T] 
which shows the comparison between the velocity distribu- 
tion of a system composed of the Milky Way-like halo plus 
the contribution at z = of one of the extra progenitors and 
the Milky Way-like halo only. 

The effect of the addition of an extra progenitor of mass be- 
tween 10 7 Mq and 10 s Mq to the merger tree of the Galactic 
halo is always negligible. The contribution is not negligi- 
ble anymore when the mass of the new merging subhalo is 



Mass Mq / Infall Time (Gyr ago) 


1.75 


2.0 


7.5 


13 


10 7 


X 


X 


X 


X 


10 s 


X 


X 


X 


X 


10 9 


/ 


X 


/ 


/ 


10 10 


/ 


/ 


/ 


/ 



Table 1. Table showing the presence (or absence) of differences 
between the velocity distribution of a system composed of the 
Milky Way-like halo plus the contribution of a single subhalo and 
the Milky Way-like halo only. Different masses are considered for 
the new progenitor, as well as different infall times. The symbol 
/ is used if there are differences in these two values of /(r, v) (in 
a range of six orders of magnitude). On the other hand, if the 
values of the two /(r, v) are identical, the symbol X is used. 



equal to 10 9 Mq. If plotted separately, without the presence 
of the "background" distribution produced by the Galactic 
halo, the contribution is present and clearly visible. Unfor- 
tunately, once we merge in the same plot these two contribu- 
tions, the one of the extra progenitor is not visible anymore 
because it is concealed by the smooth "background" distri- 
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Fig. 5a: 10 8 M Q s£ M < 10 9 M e Fig. 5b: 10 9 M < M < 10 10 M 
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Fig. 5c: 10 10 Mq^M < 10 11 M 
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Figure 5. Series of plots presenting the contribution to the Milky Way-like halo angle distribution, simulated in Fig. [4] from the satellites 
with mass in the range 10 8 Mq ^ M < 10 9 Af© (Fig. 5a), 10 9 M sg M < 10 10 M (Fig. 5b), and 10 10 M sC M < 10 11 M (Fig. 5c). 
The greyscale and the dimensions of each pixel are the same as used Fig. [4] 



bution produced by the Galactic halo. 

The only case in which the merger of an extra progenitor 
generates clear features in the angle distribution of the Milky 
Way is when the mass of the subhalo is equal to 10 10 M . 
This configuration is shown in the four snapshots of Fig. [5] 
Fig. [6^ represents the fall of the subhalo into the Galaxy ap- 
proximately 1.75 Gyr ago. Two vertical, narrow stripes are 
clearly visible. They both have an angular width of approx- 
imately 15° — 20°: the first does not have sharp contours, 
while the second is centred in the range 79° < 9 < 96° 
(—0.2 < cosf? < 0.2). The presence of these two stripes 
shows that the subhalo is beginning to get stripped and 



that streams start to form. Their width allows them to be 
potentially detectable by current directional experiments. 
At t2 the two vertical stripes have almost completely dis- 
appeared (Fig. UDd) • This behaviour is expected because the 
main body of the satellite is still coherent and it is not lo- 
cated anywhere near the Solar neighbourhood. 

At to, the satellite has performed several orbits around 
the Galactic centre (Fig. |6};) . At each orbit the phase mixing 
acts on the satellite, increasing the number of streams and 
the amount of debris into which the subhalo is disrupted. 
The streams overlap and the velocity distribution of the sub- 
halo becomes smooth and uniform. Nevertheless, there is a 
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Fig. 6a: t\ ~ 1.75 Gyr Fig. 6c: t 3 ~ 7.5 Gyr 
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Figure 6. Series of snapshots presenting the angle- velocity plot at z = of a system composed of the Milky Way-like halo and a subhalo 
of mass 10 10 M©. The four plots describe configurations in which the progenitor falls into the host halo on a radial orbit at different 
times. In Fig.[6ji it fell ~ 1.75 Gyr ago. Fig. |6p corresponds to a infall time of ~ 2.0 Gyr ago, while in Fig. [6}; (3 ~ 7.5 Gyr ago. Finally, 
in Fig. [fjjl the satellite merged ~ 13 Gyr ago. The greyscale and the dimensions of each pixel of the four snapshots are the same as Fig. 
[4] and the angle distribution is smoothed over a (5 X 5) box. 



narrow, vertical stripe at velocity ~ 430 km s 1 (v ~ 0.75), 
centred in the range 79° < 9 < 96° (-0.2 < cos6» < 0.2). 
This is due to the presence of a dense, massive stream at 
the Solar position, and it suggests that the disruption of the 
10 10 Mq satellite is not complete. 

Finally, Fig. [6ji presents the configuration in which the 
merger of the subhalo happened 13 Gyr ago. The stripes 
that are present in other plots are now completely washed 
out. The absence of stripes confirms the presence of an al- 
most homogeneous distribution. Nevertheless, the remnant 



of the central part of the satellite is not completely disrupted 
and it still shows up as an overdensity. This overdensity is 
located in the speed range 400 km s^ 1 < v < 430 km s _1 
(0.7 km s" 1 < v < 0.75 km s _1 ), with an angular width of 
about ~ 5° (-0.05 < cosf) < 0.05). 
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4 DISCUSSION 



The conclusion inferred from the analysis of the energy spec- 
trum of a system resembling the Milky Way is that the ultra- 
fine DM distribution in the Solar neighbourhood is com- 
posed of a large number of streams. The streams overlap 
and produce a smooth distribution. Our result is in agree- 
ment with the conclusion of VW, which has been reached 
using an approach complementary to this analysis. VW's 
technique must be used in tandem with a N-body simula- 
tion. This implies the need of a large computational time 
and it requires a softening length of 3.4 kpc to achieve sta- 
ble results. In our model, the gravitational force does not 
have to be artificially softened. Consequently, the phase- 
space can be explored very rapidly and accurately in a single 
timestep, without the computational overhead of numerical 
integration, and the mpc-scale resolution required by direct 
detection experiments can be reached by using a less realis- 
tic gravitational potential. The fact that these very different 
techniques give the same answer provides some confidence 
that the results are not an artifact of the approximations 
made in each case. 

When we investigate the velocity distribution /(r, v) in the 
(cos 8, v) space, the results suggest the presence of two levels 
of overdensities. They are the marks left by recent merger 
activity and we have verified that they have been mostly pro- 
duced by satellites with mass in the range lo 10-11 M . The 
probability of detecting such overdensities is however low. 
This is mainly caused by both the lack of strong features, 
and the limited angular resolution of the current generation 
of instruments. Nowadays this resolution is ~ 15° (Dujmic et 
al. 2008, Ahlen et al. 2010), one order of magnitude larger 
than that required for the detection of these overdensities 
(~ 1°). 

The analysis of the contribution from single mergers 
to the ultra-fine DM distribution of a Milky Way-like halo 
shows the presence of features when the mass of the addi- 
tional satellite is equal or larger than 10 9 Mq. The features 
generated by the merger of a satellite of mass 
are concealed by the smooth "background" distribution pro- 
duced by the Galactic halo, that dominates the (cos#,t>) 
space. Once this "background" is subtracted, it is possible to 
see the imprint left by the additional satellite. Nevertheless, 
its angular width is too small (~ 1°) to be detectable with 
the present-day technology. Moreover, the angular width of 
the features generated by the merger of a 10 10 Mq progeni- 
tor is 15° — 20°, which is accessible to ongoing experiments. 
This raises the possibility to detect these peculiar features. 
This possibility depends on the probability of finding the 
remnant of a 10 10 Mq subhalo in the Solar neighbourhood 
and that this subhalo merged into the Galactic halo in the 
recent past (< 7.5 Gyr ago). This probability is not neg- 
ligible: various works (Stewart et al. 2008, Purcell et al. 
2007, Zentner 2007) found that the mass of a DM halo, 
Mh, is provided mostly by the accretion of subhalos of mass 
~ (0.05 — 0.15)M?j,. Furthermore, it has been estimated that 
the probability of a Milky Way-like system of having experi- 
enced a merger with a subhalo of 10 11 M Q in the last 8 Gyr 
is ~ 80% (Stewart et al. 2008). 



5 SUMMARY 

In this paper we have refined the method developed in Fantin 
et al. (2008), applying it to the formation of the Milky Way 
in a cosmological context. The refinement consists of com- 
bining a merger tree, which describes the hierarchical growth 
of a Milky Way-like halo, with the model previously devel- 
oped, which simulates the interaction between an unbound 
system of particles and a larger parent halo. This allows us 
to produce a more complete and realistic treatment of the 
merger history of the Galactic halo, mapping out the sub- 
mpc-scale structure of the halo. Using such a simulation we 
can obtain new insights into the likely signature of a halo 
merger event in a small terrestrial DM detector. 
We find that the velocity distribution does not contain any 
clear feature. This is caused by the overlapping of a huge 
number of streams, which generates a uniform and smooth 
velocity distribution. This result is in agreement with recent 
works, such as Vogelsberger & White (2011) and Schneider 
ct al. (2010). Detectable signatures in the ultra-fine velocity 
distribution in the Solar neighbourhood are present only in 
the case of the recent merger into the Galactic halo of a sub- 
halo of 10 10 M Q . In this particular configuration the angular 
width of these features is larger than 10° , a resolution acces- 
sible to ongoing experiments. Finally, in the scenarios that 
we have analysed, the presence of some overdensities has 
been identified. The current generation of detectors does 
not have the angular resolution required to observe these 
features, but a future generation of detectors with resolu- 
tion of ~ 1° would resolve them, allowing the recent history 
of the Milky Way to be revealed. 
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